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Abstract 

Thermodynamic properties of the multiband superconductor MgB2 have often been described 
using a simple sum of the standard BCS expressions corresponding to a- and 7r-bands. Although, 
it is a priori not clear if this approach is working always adequately, in particular in cases of strong 
interband scattering. Here we compare the often used approach of a sum of two independent bands 
using BCS-like a-model expressions for the specific heat, entropy and free energy to the solution 
of the full Eliashberg equations. The superconducting energy gaps, the free energy, the entropy 
and the heat capacity for varying interband scattering rates are calculated within the framework of 
two-band Eliashberg theory. We obtain good agreement between the phenomenological two-band 
a-model with the Eliashberg results, which delivers for the first time the theoretical verification to 
use the a-model as a useful tool for a reliable analysis of heat capacity data. For the thermodynamic 
potential and the entropy we demonstrate that only the sum over the contributions of the two bands 
has physical meaning. 

PACS numbers: 74.25Bt, 74.70.Ad, 74.62.Dh 
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I. INTRODUCTION 



Apart from the high transition temperature of 40 K,— two-band superconductivity was 
the other unexpected phenomenon in MgB 2 which attracts increasing attention. In fact, at 
present it appears that MgB 2 is the only superconductor with substantiated theoretical and 
experimental evidence for two-band superconductivity. 

Historically, two-band superconductivity has already been investigated theoretically 
shortly after the formulation of BCS theory. Suhl, Matthias and Walker 2 suggested a model 
for transition metals considering overlapping s- and <i-bands. At the same time, Moskalenko 
proposed an extension of BCS theory for multiple bands. 3 A review of theoretical treatment 
of the critical temperature T c of multiband superconductors may be found in Ref J. 

In the early sixties there have been experimental claims for the observation of two-band 
superconductivity in some transition metals like e.g. V, Nb and Ta^ and later in oxygen 
depleted SrTi0 3 . 7 

Until now, MgB2 appears to be first system for which multi-band superconductivity has 
independently been evidenced by several experimental techniques like, for example, heat 
capacity, tunneling spectroscopy, Raman spectroscopy, penetration depth measurements, 
ARPRES, and the analysis of the critical fields. 8 The theoretical justification for two-band 
superconductivity in MgB2 has been given from electronic structure calculations.^^ These 
find that the Fermi surface contains two quasi-cylindrical sheets corresponding to nearly 
two-dimensional a-bands. A three dimensional network of the Fermi surface is attributed 
to the 7r-bands. It has been demonstrated that the optical bond stretching E 2g phonons 
couple strongly to the holes at the top of a-bands while the three-dimensional 7r-electrons 
couple only weakly to the phonons. The different coupling strengths of the a- and 7r-bands 
lead to superconducting gaps different in character and size . n 4 2 i 13 i 14 i 15 Using linear response 
theory it is possible to calculate the electron-phonon coupling (Eliashberg functions) from 
first principles. The superconducting gaps obtained from Eliashberg theory are in very good 
agreement with the experiments . l16 

Interband scattering from impurities will complicate this picture because interband scat- 
tering leads to a decrease of T c and finally to a single order parameter . 16 i 17 i 18 i 19 Interband 
scattering is weak in MgB2,— but this is not necessari ly the case in samples in which Me 



has been replaced by Al or B by C (Ref|2j2j2j2Jj2 



26.27.28.29.30.31.32.33.34.35.36.3 
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or which have been exposed to neutron irradiation.— Such samples exhibit considerably 
reduced T c , while the two gaps persist even at very low critical temperatures. Recently it 
was shown that the T c reduction in MgB 2 due to Al or C doping can be explained mainly 
as due to a simple effect of band filling. 39 ' 40 A similar observation has been made using 
a phenomenological weak-coupling approach.— Further, the doping independent 7r-gap in 
C-doped MgB2 can be understood as due to a compensation of band filling and interband 
scattering effects. 

Thermodynamic properties of anisotropic superconductors in the weak coupling regime 
were extensively studied in the past. In the case of weak anisotropy the BCS model was 
extended by Pokrovsky. 42 It was shown that the specific heat jump at T c is reduced as 
compared to the isotropic case. For two-band weakly coupled superconductors the specific 
heat was calculated by several au thors 43 ' 44 ' 45 ' 46 ' 47 ' 48 ' 49 ' 50 ' 51 ' 52,53 ' 54 (for a recent review sec also 



Ref 



55T ). The main prediction is that at T c the relative jump in the electronic specific heat, 



{Csc — Cn)/Cn, is reduced as compared to the universal BCS value of 1.43. On the other 
hand, for an isotropic strongly coupled superconductor the relative specific heat jump is 
larger than 1.43 (see e.g. the review in Ref. 56). The combined effect of strong coupling and 
multiband anisotropy on the specific heat was studied earlier by the present authors,— where 
the results of the first principles calculations of the electron-phonon interaction in MgB2 were 
used but the effect of interband impurity scattering was not considered. Recently strong- 
coupling corrections were taken into account in the so-called two-square-well approximation 
(separable model) , ?7iWi?9 where the effect of interband scattering on some thermodynamic 
functions was studied . 5Sl59 

In the present work we formulate a generalized description of the thermodynamics of 
multiband superconductors taking into account impurity scattering (magnetic and nonmag- 
netic) in the framework of two-band Eliashberg theory. The results are applied to MgB2 
using the first principles band-structure results for the electronic spectra and electron-phonon 
interaction^ 2 - by extending our preceding approach. 16 The superconducting energy gaps, the 
free energy, the entropy and the heat capacity for varying nonmagnetic interband scattering 
rates are calculated within the framework of two-band Eliashberg theory. It will be shown 
that the expression for the thermodynamic potential on the extremal trajectory correspond- 
ing to solutions of the Eliashberg equations has the form of the sum of contributions of a- and 
7r-bands, but that only the total thermodynamic potential (the sum of both contributions) 
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has physical meaning. 

In a second step, we perform a comparison of the phenomenological two-band a-model 
with the Eliashberg results and apply a fit program developed for the a-model to extract the 
gaps and the Sommerfeld constants from the Eliashberg results. Good agreement of the two 
band a-model with the Eliashberg data is found for the temperature dependence of total 
heat capacity, the entropy and the free energy and the gaps. There are, however, distinct 
deviations in the partial contributions to the individual quantities and the Sommerfeld 
constants obtained from the fits. We conclude that the phenomenological a-model approach 
can be taken as a handy tool to analyze e. g. experimental heat capacity data and the gaps 
to a satisfying accuracy, however that care must be taken for other quantities. 

The paper is organized as follows: In Section II the introduction to the formalism and the 
method of solution is given, in Section III numerical results for the densities of states (DOS) 
and various thermodynamic quantities as a function of interband impurity scattering rate 
are discussed, in Section IV the comparison of the two-band a-model with the Eliashberg 
results is performed. In the Appendix a general expression for the thermodynamic potential 
of a multiband superconductor with nonmagnetic impurities is derived. 

II. FREE ENERGY AND ELIASHBERG EQUATIONS 

A general expression for the difference of free energies A'f2 = Qn — ^s in the normal (N) 
and superconducting (S) state for a system with electron-phonon interaction and multiple 
bands can be obtained in two ways: one has been derived by a straightforward integration 
over the electron-phonon interaction constants by Golubov et a/..— The derivation of the 
expression for the thermodynamic potential for the case of nonmagnetic as well magnetic 
impurities is presented in Appendix A. In terms of Matsubara frequencies the f2-potential 
can be written as 



where Qe is the f2-potential of the noninteracting electrons, and f2p h is the fi-potential of 
the noninteracting phonons. 
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For the difference of the f2-potentials in the normal and the superconducting state one 
obtains 



AQ = Q N -tt s = -7rJ2 N i(°) T Yl 



(1) 



n=—w, 
, 2 



7N 



D- 



^l[{Zf n ) -l] + 2$f n , JiZUZi - 1) + $ 



Z is the renormalization factor (which is unity in the weak coupling limit) and $ is the order 
parameter which is connected to the energy gap via Aj(u; n ) = $j(k n ) / Zi(cu n ) = § in /Z in . Z N 
and Z correspond to the normal state (A = 0) and the superconducting state, respectively. 
The summations in Eq. Q are carried out over the fermionic Matsubara (temperature) fre- 
quencies u n = irT(2n — 1) as well as over the band index % = a, tt. iVj(O) are the partial 
electronic DOS's for the a- and 7r-bands at the Fermi level. 

Another way to express the free energy was suggested by Carbotte,— who proposed 
that it is possible to find a functional, which minimization with respect to Z and $ gives 
the Eliashberg relations for superconductors with strong electron-boson interaction. For a 
multiband system the corresponding functional is given by 



+ 



in J 
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\Au n -u m ) - n*Au c ) , (2) 



■ jm 



where 



du 2 a%Au)Fij(u) 



+ (Tij ± r f 7) 5 n , m /27rT 



represents the electron-phonon interaction together with nonmagnetic and magnetic 
T™j impurity scattering terms. The Coulomb pseudopotential fi*Au c ) is determined at a 
frequency uo c which has to be chosen much larger than the maximal phonon frequency. 
Minimization of Eq. [21 provides the Eliashberg equations on the imaginary (Matsubara) axis 
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Z ni 0J n = U n + 1lT^2 >Hj{u n - 0J rn ) Z Jrn UJ m ^ 



m,j y (Z jm UJ n ) 2 + $ 2 m 

l^ml<Nc| - 



- ni 



jm 



(Z jm oj m ) +$| r 



Both functionals AQ and AjF give the same result on the extremal trajectory which 
corresponds to solutions of the Eliashberg equations Eq. (jHJ) • In the following we will use AQ 
for the calculations of thermodynamic quantities. 

For the calculations e.g., of the densities of states Ni(u) in the superconducting state, 
we need to know the renormalization factors Zj and the order parameters $j along the real 
frequency axis. The corresponding analytical continuation of Eq. (jH)) substituting ioj n ==>- 
oj + i5 gives 



ZAoj)oj = oj 



POO 

doj' / dvc^-^Fijiy^Iiyj + i5, u, u/)Re 
„■ Jo 



Zj(u')L0' 



~7 J( > J(Z 3 (oj')oj') 2 - <£> 2 (u/) 
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/oo poo 
du'S^ / dvo%Av)Fij{v)I{u + i8,v,uj')R£- 



4/2 J>*-(u; c ) jT ' dw'tanh Re 



;^VV') 2 -^V) 



+ 



where 



E - ^) y — ( 5 ) 



i (a; + zd, z/, oj ) = + 



uj + ib — v — a;' c<j + i5 + z/ — oj' 
and n(z/) and /(a/) are Bose and Fermi distribution functions, respectively. Please note that 
the nondiagonal elements of all functions o^j{y)Fij{y), fj,*JoJ c ), T^and T 7 ^ have to satisfy 
the requirement of the detailed balance principle 

N a (0)T a „ = N n (0)T Ta , 
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where N a (0) and A^-(O) are the normal state bare electronic DOS's in the a- and the 7r-bands, 
respectively. 



As has been shown in Ref. 



581 ando^. for a separable interaction Xij(u> n — u m ) = Xij6(Q — 
\u) n \)9(Q — \uj m \) the standard weak coupling expressions can be obtained which correspond 
to the two-band BCS results. 

In this paper we will use the following representation of the Eliashberg equations which 
is better suited for numerical solution by iterations^ 
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(6) 
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\2T 



Re 



(7) 



where $j(u;) and = Zj(u;)u; are the renormalized gap function and the renormalized 

frequency respectively, denotes the impurity scattering rate within the Born approxi- 
mation. The real and the imaginary parts of the Eliashberg functions and u)i(u) are 
connected by the Kramers-Kronig relations. Hence, they have the same Fourier images. This 
yields a procedure for a rapid solution. The convolution type integrals (Eqs. EHZ|) should 
be calculated by the Fast Fourier Transform (FFT) algorithm. The inverse complex Fourier 
transformations of the results obtained give complex values of $j(u/) and Ui(tj). 



III. INTERBAND SCATTERING 



Intraband scattering from nonmagnetic impurities does not affect T c and the supercon- 
ducting densities of states Ni(u), as well as the thermodynamic potentials. However, inter- 
band scattering is expected to modify T c and NAuj) strongly. In the weak coupling regime 
this effect has been demonstrated in Refsl^ Ml^jlsl In the following we will calculate T c , the 
gap functions and the superconducting DOS by solving the nonlinear equations (Eq. |3]-[7J) for 
various values of the interband nonmagnetic scattering rates T aix and f ro . For convenience, 
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we define an interband scattering parameter T in the following way T U7T = TN n (0)/Ntot(0), 



r w = TN a (0)/N tot {0). 



A. Gap Functions and the Density of States. 

As shown in Fig. ^ T c gradually decreases with increasing T and saturates at a value 
corresponding to that expected for isotropic coupling. The initial decrease of T c with F 
amounts to T C (T) - T c (0) = 0.10(1) K/cnT 1 ■ T, 




100 200 300 400 500 

r (cm" 1 ) 



FIG. 1: The dependence of T c on the interband scattering parameter T. 



The DOS in the superconducting state, Ni(u), is given by the expression 

ZAu) u 



N t (u) = A,(0)Re- 



(8) 



where Aj(0) is the DOS in the normal state at the Fermi level of the corresponding energy 
band. Here 3>i(o>) = Ai(u)Zi(u), where Aj(o;) and Z{{uS) are complex pair potentials and 
renormalization functions. Figs. El and El display the Z a ^{oj) and A tT 7r (u;) as obtained from 
using spectral functions calculated from first-principles for the effective two-band model in 
MgB 2 .i& 

The results demonstrate the self-energy effects arising due to the sizeable electron-phonon 
interaction in MgB 2 . The real parts A CTj7r (u;) and Z a ^{io) strongly depend on uj when uj be- 
comes comparable to the characteristic phonon frequencies. The imaginary parts appearing 
at these energies indicate the decay of quasi-particles due to this strong interaction. Fur- 
thermore, the effects of impurity scattering are also visible as additional structure at low 
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FIG. 2: The energy dependencies of the renormalization functions in a two-band model at T/T c = 
0.1 in the clean limit T = and for r = 50cm -1 . 

energies comparable to the scattering rate T. This structure is particularly strong in the 
real and imaginary parts of Zi(u). The latter can be seen from the last term in Eq. 0] The 
impurity contribution to ImZj(o;) is proportional to riVj-(w), where i,j belong to different 
bands. 

Fig. |3] shows the densities of states for different magnitudes of the interband scattering 
rate T at low-temperature (T/T c = 0.1). In the clean limit, the two bands show two 
different excitation gaps. In accordance with earlier calculations , 17118 the interband impurity 
scattering mixes the pairs in the two bands, so that the states appear in the cx-band at the 
energy range of the 7r-band gap. These states are gradually filled in with increasing scattering 
rate. At the same time the minimal 7r-band gap in the DOS raises due to increased mixing to 
the cr-band with stronger electron-phonon coupling. Thus the decrease in T c is accompanied 
by an increase of the minimal gap in the excitation spectrum as has been observed by 



Gonnelli < t <il. in kefl^ and theoretically supported by some of us in Rof. 



39 



Fig. El shows the evolution of the superconducting DOS with temperature for a fixed 
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FIG. 3: The energy dependencies of the gap functions obtained by using a strong coupling two-band 
Eliashberg model at T/T c = 0.1 in the clean limit and with r = 50cm~ 1 . 

values of the interband scattering rate T = 10 cm -1 . One can see that at finite temperature 
the densities of states in both bands become gapless: In addition to the states at the 
energy range between the 7r-band and the cr-band gap, states appear down to lowest energies 
due to thermal phonons. Such gapless behavior is most pronounced close to T c . In the 
isotropic single-band superconductor, this thermal effect in the strong-coupling regime was 
demonstrated earlier in Reflfill Note, that the shape and temperature dependence of the 
superconducting DOS are very different compared to the sum of two BCS-like densities 
of states. This is particularly pronounced for the cr-band. Therefore, one would expect a 
non-BCS temperature behavior in the thermodynamical functions. 

B. Thermodynamic Functions 

For the numerical calculations for the free energy difference we have used Eq. (0) with 
parameters described in the previous Section. The entropy difference between the normal 
and superconducting states is determined by the first derivative of the free energy difference 
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FIG. 4: The evolution of the low-temperature densities of states with the interband scattering rate 
r in a strong coupling two-band Eliashberg model at T/T c = 0.1. The upper panel a) shows the 
superconducting DOS for the u-band, the lower b) the same for the 7r-band. 

with respect to temperature 

AS{T) = dAtt{T)/dT, 
and the specific heat difference by the second derivative with respect to temperature 

AC(T) = Td 2 AQ{T)/dT 2 . 

Here we note that taking derivatives from numerically calculated data (as well as from 
experimental ones) is often a mathematically ill-defined or numerically unstable procedure. 
Therefore, we used three different schemes to interpolate the numerical data: a) a Chebyshev 
scheme to interpolate the free energy calculated at non-equidistant points Tj = cos(^), (j = 
0,1, n; where n = T c j AT is the number of points) and constructing a corresponding matrix 
nxn operator,— b) a polynomial approximation which works well for large temperatures 
where the densities of states are smooth functions without square-root singularities, and 
c) an interpolation of the free energy differences by a series of Bessel functions Ki(nA/T) 



similar to the weak-coupling BCS approximation (see e.g., Ref.63) . The latter captures 
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FIG. 5: The evolution of the superconducting densities of states with temperature at fixed interband 
scattering rate T = 10 cm -1 in a strong coupling two-band Eliashberg model. 

the superconducting square-root features and works well at low temperatures. The data 
presented below were chosen such that all three approaches gave similar results. 

The specific heat jumps AC(T C ) at T = T c were determined separately by the calculation 
of the coefficient (3 = T c AC(T c )/2 in the term Afi(T -> T c ) = (3(1 - T/T c ) 2 . 

First, we consider the case without interband impurity scattering ('clean' case T = 0). 
The expression for AQ(T) (see Eq. EJ) consists of two terms containing iVj(O), the renormal- 
ization factor Zj, and the order parameter $, (or the energy gap Aj = $j/Zj) for each band 
separately. These terms reflect the partial contributions of each band to the total free energy 
(cf. upper panel of Fig. [BJ). One sees that the 7r-band gives a negative contribution to the 
free energy over the full temperature range. This surprising observation reflects the fact that 
creation of the superconducting state in the 7r-band and, since the superconducting state 
in the 7r-band is induced by the occurrence of superconductivity in the strongly interacting 
cr-band, costs energy. 

The coupling due to the nonzero off-diagonal elements of the electron-phonon interaction 
is the reason for the same critical temperature T c and the induced superconducting order 
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FIG. 6: Partial and total contributions to the free energy (AO), the entropy (AS) and the heat 
capacity AC/T in the clean case. 

parameter in the 7r-band. But the analogy to two-independent contributions to the free 
energy is not fully applicable. The partial Aflj's near T c behave as 0(T C — T) instead of 
AQ tot oc 0((T C — T) 2 ), according to the requirements of a mean field theory. In the middle 
panel of Fig. |U] one can further see that the entropies AS , (Ti7r (T c ) have finite values, whereas 
only AS to t{T c ) = 0, as required by the third law of thermodynamics. According to this one 
has to consider only the total thermodynamic functions as physical ones. 

Effects of interband impurity scattering on thermodynamic functions are shown in Fig. 
The reduced specific heat jumps AC(T C ) /T c at T = T c grows monotonically with the increase 
of r from 3.24 mJ mol _1 K~ 1 in the clean case^ to 4.1 mJ mol _1 K _1 for T =1000 cm -1 . These 
values correspond to ratios AC(T c )/'~f(T c )T c « 0.98 (clean case) and 1.3, which are smaller 
then the corresponding BCS value of 1.43 in a single band model. At low temperatures the 
ratio AC(T)/T saturates to the value \im.T-,o(Csc(T) /T — Cn(T)/T) = — 7a?(T = 0) = 
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FIG. 7: Total contributions to the free energy (Af2), the entropy (A 5) and the heat capacity 
AC / T for various impurity scattering rates T. 

—7o(l + A a „) = —3.24 mJ mol _1 K -2 which is determined by the bare (band) electronic 
specific heat capacity 70 = 27r 2 fc 2 3 (iV o .(0) + N n (0))/3, and the average coupling constant 

A N a (0)(X aa + X an ) + N n (0)(X wlT + A^r) 
N a (0) + N v (p) 

and does not depend on impurities. 

IV. ELIASHBERG VERSUS TWO-BAND a-MODEL 

Since the two-band a-model has been widely used to analyze the experimental heat 
capacity data of MgB2 it was interesting to see to which extend a two-gap a-model can 
reproduce the Eliashberg results and, if so, how the corresponding parameters compare with 
those identified from the Eliashberg calculations. 
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The a-model originally introduced by Padamsee et al, in close analogy to the BCS theory, 
assumes a BCS temperature dependence of the superconducting gap. The magnitude of 
superconducting gap at T = is introduced as an adjustable parameter a (from which the 
a-model received its name).— The parameter a is defined according to 

A(T) = (a/«Bcs)A BC s (9) 

with «bcs — Abcs(O) / ksTc = 1.764 being the weak-coupling value of the gap ratio. 
Within the scope of the a-model the free energy F$ in the superconducting state can be 
written as 

oo 

fis(T) = 2 N(0) J de {2 k B T ln[l - f(E)} - \ - } (10) 



with E = \/e 2 + A 2 and A^(0) being the electron and phonon renormalized band-structure 
electronic density of states at the Fermi energy. 

We subtract the normal state contribution to the free energy which corresponds to A = 
and introduce the dimensionless parameters t = T /T c , 5(t) = A(T) = Ao and x = e/ ksT c 
and arrive at 



00 

Am = 2N( ) (far.)' ftoht H^£$] - 5 - a2d(i)V(!, ' a) } 

J I 1 - f\ x ) 2 y y J 



(11) 

where y = a/x 2 + a 2 5(t) 2 , f(y,a) = l/[l + exp(y/*)], and f(x) = f(y,a = 0) = 
l/[l+exp(x/t)] 

The electronic entropy in the superconducting state is obtained from the first derivative 
of the free energy with respect to temperature and can be written as 

S , el (t)/ 7 T c = -(3/7r 2 ) J dx{f(y,a)lnf(y,a) + (1 - f(y,a))ln(l - f(y,a))} (12) 



wherein the normal-state electronic specific heat capacity ('Sommerfeld-term') is given 
by 7 =|iV(0)7r 2 k B . 
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The electronic heat capacity C e \ is calculated from Eq. (|12|) according to 

C d /yT c = t(d/dt)S el / 1 T c (13) 

To compare the Eliashberg results with the two-gap a-model approximation we have 
developed least-squares refinement codes to fit the entropy (Eq. H2J) and the heat capacity 
(Eq. IT3*|) with an a-model which linearly superposes the contributions from the a and the ir 
electronic system. For the temperature dependence of the reduced gap S(T) = A(T) / A(0) 
we adopted the tabulated values provided by Miihlschlegel.— For the analytical calculations 
we used a polynomial fit of the these data. 

The heat capacity was calculated from Eq.H5l using an appropriate numerical difference 
quotient as approximation for the derivative with respect to t. Integrations in Eq. (jl2j) 
were performed numerically with a Gaussian quadrature scheme with a cut-off for x> 100. 
Examples for the fits of the entropy and the heat capacity for T = and T = 1000 cm -1 
corresponding to T C = 39.4K and T C = 26.5K, respectively, are displayed in Fig.[HJ 

To fit the data we varied aij viz. the energy gaps Aj(T) and the Sommerfeld constants 7, 
(i = a, 7r) and a single critical temperature T c . All results are compiled in TableQJ Fig.lTUl 
displays the fitted gaps versus T c viz. the interband scattering parameter T. 

Attempts to fit also the total energy were less successful and provided results inconsistent 
with the results of the entropy and heat capacity fits. In these fits we observed a tendency to 
converge to essentially a single-band model with an averaged gap somewhat above the weak 
coupling BCS result of A(0) / &bT c = 1.76. Calculations of the free energies with the on and 
7i parameters obtained from the fits to the entropies and heat capacities reproduced the 
Eliashberg free energies equally well as obtained from the fits of Eq.EB A closer inspection 
revealed that characteristic differences for various T are only visible at small temperatures 
(< T c /5) where the free energy levels to saturation. The fits apparently are not sensitive 
enough to catch these slight deviations at a satisfying level. 

Table|l] compiles the parameters obtained from the fits of the entropies and the heat 
capacities. The results are largely independent of whether they are obtained from fits of the 
entropy or the heat capacity. In general, gaps obtained by fitting of the entropy are closer 
to the Eliashberg (Matsubara) gaps. Fits in the clean case (r — ► 0) readily converged with 
the parameters listed in Table |1J For large values of T, convergence of the fits of the heat 
capacities with the two-band model were less stable and fits with a single-band model in 
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FIG. 8: (solid lines) Fits of the Eliashberg entropies and the heat capacities for T = (T C =39.4K) 
and r = 1000 cm -1 (T c =26.5 K) with the a-model (Eq. El an d Eq. 1131 For both entropy data sets 
(r = cm -1 (o) and 1000 cm -1 (□) and the heat capacity data set with T=0 cm" 1 (o) a fit with 
the two-band a-model was carried out. The fit of the heat capacity data for T= 1000 cm -1 (□) 
was performed with a single-band a-model alone. 

some cases proved to be more conclusive. 

In Fig. Elwe show the total and the partial contributions to the free energy, entropy and 
the heat capacity calculated according to the a-model using the fitted parameters given in 
Table H] for T = (T c = 39.4 K). The a-model describes the total heat capacity rather well. 
There are subtle differences in the partial 7r and a contributions below T « 10 K. These 
difference are also reflected in the fitted ratios of the Sommerfeld constants (ja/l^&t which 
deviates markedly from the ratio of the phonon renormalized Sommerfeld terms used for the 
Eliashberg calculations, (7^/7^) = N^/N^O) ((1 + A CT(T + A CT7r )/(l + A™ + A^)) « 1. 

Naturally, since within the scope of the a-model all partial contributions are positive 
definite, the negative 7r partial free energy and the sign change of the ir partial entropy 
(compare to Fig|UJ) cannot be reproduced. 

Finally, Fig. El shows the superconducting gaps as obtained from the fits in comparison 
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T c (K) 



r (cm" 1 ) source Af iash (meV) A a (meV) Af iash (meV) A^ (meV) 7^/7^ 



39.4 



38.1 



36.9 



35.5 



34.0 



32.5 



30.1 



26.6 



10 



23.8 



40 



64.4 



100 



196 



1000 



s 
c 
s 
c 
s 
c 
s 
c 
s 
c 
s 
c 
s 
c 
s 
c 



7.04 



6.51 



6.12 



5.77 



5.47 



5.95 



4.80 



4.23 



7.48 
7.16 
6.97 
7.00 
6.71 
6.39 
6.34 
6.30 
5.68 
4.57 
5.45 
4.51 
5.02 
4.49 
4.44 
4.30 



2.67 



2.92 



3.14 



3.36 



3.56 

1.38 
3.57 



3.97 



4.09 



2.88 
2.51 
3.12 
3.15 
3.40 
3.05 
3.58 
3.58 
3.63 



3.92 



3.59 



0.72 
0.74 
0.63 
0.62 
0.64 
0.67 
0.53 
0.52 
0.66 

0.60 

0.80 

2.7 



TABLE I: Critical temperatures T c , Aj and Sommerfeld parameters 7^ as obtained from least- 
squares fit of the polynom interpolated Eliashberg entropies (S) and the heat capacities (C) by 
a two-band a-model. T is the interband scattering parameter. The sum 7 = j a + j n of the 
fitted Sommerfeld terms was found to be constant within 3%. For comparison the Eliashberg 
(Matsubara) gaps (solutions of Eq. as displayed in Fig. E3are listed. 

with the Eliashberg calculations. The agreement is fairly good for higher T c . Deviations are 
seen for T = 33 K for the gaps gained from the fits of the heat capacities, while the gaps 
received from the fits of the entropy rather well follow the Matsubara calculations and the 
merging point of both gaps at the weak coupling value is also well reproduced. 
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FIG. 9: Partial and total contributions to the free energy (Af2), the entropy (AS) and the heat 
capacity AC/T in the clean case calculated with the two-band a-model. 

V. CONCLUSION 

In summary, using the Eliashberg approach, we have studied the behavior of the su- 
perconducting density of states, energy gaps, free energy, entropy and specific heat in a 
strongly-coupled two-band superconductor with interband impurity scattering. We have 
demonstrated strong modifications of the densities of states by interband scattering and 
have shown how thermal effects modify these results. We have calculated the temperature 
dependencies of the free energy, the entropy and the specific heat and the specific jump 
at T c as a function of interband scattering rates and performed a detailed comparison of 
the phenomenological two-band a-model with the Eliashberg results. We have shown that 
despite strong modifications of the DOS by interband scattering, the a-model approach is 
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FIG. 10: Superconducting gaps A(0) CT and A(0) 7r as obtained from the fits of the Eliashberg 
entropies (o) and the heat capacities (□) with the two-band a-model (Eq. 1121 and H^)l . Gap data 
gained from fits of the heat capacities for T c < 35K were obtained from a fit with a single one-band 
a-model. The (red) solid line shows the weak coupling result Abcs = 1-76 T c . The dashed (blue) 
lines represents the results of the Matsubara calculations of A a (uj n = ttT) and A n (u n = ttT). 

sufficiently accurate and can - as a first approximation - be used to extract gap values from 
experimental heat capacity or entropy data. 

Interband scattering alone, however, is not sufficient to model the decrease of T c observed 
for Al and C doped samples Mgi_ x Al x (Bi_j / C J/ )2. As demonstrated recently, the decrease of 
T c can rather be understood in terms of a band filling effect due to the electron doping by Al 
and C and a concomitant scaling of the electron-phonon coupling constant A by the variation 
of the density of states as a function of electron doping.— Compensation of band filling and 
interband scattering effects shifts the merging point of the a and it gaps to higher doping 
concentrations and lower T c 's than expected based on interband scattering considerations 
alone. Only the combination of interband scattering with band filling effects allowed us 
to model the nearly constant 7r gap and the decreased critical temperature and increased 
doping concentrations at which the a and n gaps finally merge. 
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APPENDIX 

In order to calculate the thermodynamic potential Q for a superconductor with strong 
electron-phonon coupling and nonmagnetic impurities, we use a general expression for the 
electron Matsubara 2x2 matrix Green function 

A ( \ iZ jn"nT Q + Q(p)r 3 + $(p)f 2 

where £j(p) is the bare spectrum (p = {j,p,u n }, with band index j and momentum p). 
Pauli matrices f correspond to Nambu space. This Green function obeys the Eliashberg 
equations, which allows us to express the potential Q directly through it as will be shown 
below. 

The thermodynamic potential Q can always be expressed by the electron Green function 
by means of integration over the electron charge. Using e.g. Eq.(16.9) from RefQI one 
obtains 

1 . r a / A A 

G (o)(p) G(x,p) - G (0 )(p) 



fi = ^ 0) + < } +TEjf tr 



o 

^ ( e° ) +< ) +T£jftr 





A A 

E(x,p)G(x,p) 



(15) 



A A 

where x is a dimensionless factor. G(x,p) and are the exact electron Green function 

A 

and the self-energy for the case when the electron charge has the value xe. G(o)(p) is the 
Green function for zero coupling constant. 

The electron-phonon contribution can be expressed in terms of electron and phonon Green 
functions by means of Eliashberg equation: 

&(x,t M ,r) = T SK,p,j>^™^ = x 2 f 3 G(x, t M , r)f 3 D(t M , r), (16) 

<^n,P,j=T,(T 

where D(tM,r) = fl ,2 -D(o)(^M, r) + D imp (tM,r)- -D(o)(^M,r) is the phonon Green function 
expressed in coordinate representation,— 1/T < tM < 1/T is the Matsubara time. Here 
we suppose that the phonon Green function is independent of the coupling constant in the 
adiabatic approximation. This is the usual approximation, which is related to the fact that 
the electron-phonon Hamiltonian contains the phonon spectrum already renormalized due 
to the electron-phonon interaction and one should not take this renormalization into account 
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once more. The second term corresponding to impurity scattering is considered in the Born 
approximation, where D imp {tM,r) a e 2 . Below we follow Refl66l. Making use of Eq. (JTfiJl 
we can derive the simple identity: 



tr | E(x)G(x) - ~jj^G{x) 



, xdG(x)*. . 



which allows us to rewrite Eq. (|T5|) in the coordinate representation (r = {t/^,r}) 



= fi ( e 0) + SI™ + J dr J — tr £(x)G(x) - _-J^G(a;) + - ' 



2 <ix 



G(x) 



n? + ng - / */ | tr f^sw + / * / f tr f^GW 



\ / 

we can now perform the integration over x exactly and find after Fourier transformation the 
required expression: 

n = Q^+Q^+Tj2^detG^- j dr (in det G (0) ) - | £ tr (g (0) £ - f) . (17) 



This expression is valid provided that the Green functions G satisfy the Eliashberg equations. 

For the difference of free energies in the S and N-states we have after the Fourier trans- 
formation 



AQ = T ln 



A 

det Gj 
det Gf 



T 

2 



where is given by 



(z>„) 2 + e?(p) 

Eq J18I is the sum of Green functions in different bands. Finally AQ can be expressed as 



(Z/ B o; n ) 2 + g(p) + ^ n 

(^.) 2 + e 2 ( P ) 



+ 



+ r EE 



(^n) 2 + tf(p) + *? B (^>n) 2 + ^ 2 (P) 
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After integration with respect to the momentum we obtain the expression 



l^cl 



AQ = -27^53^(0) V(^n) 2 + ^ n -|^>n 

l w c| 

+7TT 53^(0) 



+ 



jrv~n 



J jn UJ n) 



I 7^ i 



This expression does not contain any impurities directly. The effect of intraband impurities 
cancels from Eqs. (3-7). Also Z s , Z N , and $ do not depend on intraband scattering, however 
these functions are dependent on interband impurity scattering. 
The final answer is the expression given in Eq. (1) 



OJr, 



j=a,7v n=—uj c 



K| + y/(zf n u n y + $ 2 n 
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